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Perturbative and Non-Perturbative Analysis of the 3’rd Order Zero Modes 
in the Kraichnan Model for Turbulent Advection 

Omri Gat, Victor S. L’vov and Itamar Procaccia 
Department of Chemical Physics, The Weizmann Institute of Science, Rehovot 76100, Israel 

The anomalous scaling behavior of the n-th order correlation functions Tn of the Kraichnan model 
of turbulent passive scalar advection is believed to be dominated by the homogeneous solutions (zero- 
modes) of the Kraichnan equation BnTn = 0. In this paper we present an extensive analysis of the 
simplest (non-trivial) case of n = 3 in the isotropic sector. The main parameter of the model, 
denoted as C,h, characterizes the eddy diffusivity and can take values in the interval 0 < ((ji < 2. 

After choosing appropriate variables we can present computer-assisted non-perturbative calculations 
of the zero modes in a projective two dimensional circle. In this presentation it is also very easy to 
perform perturbative calculations of the scaling exponent C,s, of the zero modes in the limit C,h 0, 
and we display quantitative agreement with the non-perturbative calculations in this limit. Another 
interesting limit is C,h —> 2. This second limit is singular, and calls for a study of a boundary layer 
using techniques of singular perturba tion theory. Our analysis of this limit shows that the scaling 
exponent vanishes like \JQij log Qi. In this limit as well, perturbative calculations are consistent 
with the non-perturbative calculations. 


I. INTRODUCTION 

The Kraichnan model of turbulent passive scalar ad¬ 
vection ||^ pertains to a field T{r,t) which satisfies the 
equation of motion 

+ u{r, t) ■ WT{r, t) = nV'^T{r, t) + ^{r, t). (1) 

Here is a Gaussian white random force, k is the 

molecular diffusivity and the driving field u(r,t) is cho¬ 
sen to have Gaussian statistics, and to be “rapidly vary¬ 
ing” in the sense that its time correlation function is pro¬ 
portional to 6 {t). The statistical quantities that one is 
interested in are the many point correlation functions 

Y2„(ri,r2, ...,r2„) = {{T{ri,t)T{r 2 ,t )... T(r2„, t))), (2) 

where double pointed brackets denote an ensemble av¬ 
erage with respect to the (stationary) statistics of the 
forcing and the statistics of the velocity field. Assuming 
that these correlation functions are scale invariant one is 
interested in the scaling (or homogeneity) exponent C, 2 n 
of which is defined by 

^2n(Ari, Ar2 ... Ar2„) = A‘^^"Y2„(ri, r2 . ..r^n) ■ (3) 

One expects such a scale invariant solution to exist in 
the inertial range, i.e. when all the separations sat¬ 
isfy rj Tij <C L where 77 and L are the inner and outer 
scales respectively. It is known 0 that for T 2 such a so¬ 
lution exists with ^2 = — Ch, where (^h is the exponent 

of the eddy-diffusivity, see Eq. #■ 

The Kraichnan model is unique in the field of turbu¬ 
lence in that it allows the derivation Q of an exact dif¬ 
ferential equation for this correlation function, 


[-«^X^«+^ 2„]^2n(ri,r2,...,r2„) =RHS. (4) 

OC 

The operator B 2 n = J 2 T>i 3 Ba/s, and Ba /3 are defined by 

Bafs = S(r«,r/3) = h,j{ra - r0)d‘^fdra^idr^j , (5) 

where the “eddy-diffusivity” tensor hij{R) is given by 

YAR) = h{R)[{Ch +d- 1)A, - ChR^RJ|R^] , (6) 

and h{R) = 7J(i?/£)^'*. Here C is some characteristic 
outer scale of the driving velocity held. The parameter 
that can be varied in this model is the scaling exponent 
Qh] it characterizes the R dependence of hij{R) and it 
can take values in the interval [0,2]. Finally, the RHS 
in Eq. (^) is known explicitly, but is not needed here. 
The reason is that it was argued that the solutions of 
this equation for n > 1 are dominated by the homoge¬ 
neous solutions (“zero-modes”), in the sense that deep 
in the inertial interval the inhomogeneous solutions are 
negligible compared to the homogeneous one. Also, it 
was claimed that in the inertial interval one can neglect 
the Laplacian operators in Eq. (^), and remain with the 
simpler homogeneous equations B 2 n^ 2 n = 0. 

Exact solutions of these homogeneous equations are 
not easy; even in the simplest case of n = 2 the function 
Y4 depends on six independent variables (for dimensions 
d > 2 ), and one faces a formidable analytic difficulty for 
exact solutions. Accordingly, several groups considered 
perturbative solutions in some small parameter, like A 
Q or the inverse dimensionality 1/d [Q. The rationale 
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for this approach is that at = 0 and d oo one ex¬ 
pects “simple scaling” with (^2n = ^(^2- The exponent (^4, 
and later also the set (2n, were computed as a function 
of near these simple scaling limits. The other limit 
of C/i 2 invites perturbation analysis as well, since 
one expects that at = 2 all the scaling exponents C2n 
would vanish. Such a perturbation theory turned out to 
be elusive. 

Recently we reported Q that it is possible to solve 
exactly, eigenfunctions included, the homogeneous equa¬ 
tion satisfied by the 3 ’rd order correlation function 
jF3(ri,r2,ra). Note that in Kraichnan’s model all the 
odd-order correlation functions J-2n+i are zero because 
of symmetry under the transformation T —> —T. This 
symmetry disappears for example if the random force 
is not Gaussian (but J-correlated in time), and in 
particular if it has a non-zero third order correlation 

T>3(T’i,r2,r3) = Jdtidt 2 {^iri,ti)^{r 2 ,t 2 )C{r 3 , 0 )). ( 7 ) 

With such a forcing the third order correlator is non-zero, 
and it satisfies the equation 

^3d-3{ri,r2,r3) = Vs , B3 = B12 + B13 + B23 ■ (8) 

This equation pertains to the inertial interval and ac¬ 
cordingly we neglected the Laplacian operators. We also 
denoted V3 = limr^^^o T’3(ri, r2, r3). The solution of 
this equation is a sum of inhomogeneous and homoge¬ 
neous contributions, and below we study the latter. We 
will focus on scale invariant homogeneous solutions which 
satisfy iF3(Ari, Ar2, Ar3) = A‘>^iF3(ri, r2, r3). We refer to 
these as the “zero modes in the scale invariant sector”. 
We note that the scaling exponent of the inhomogeneous 
scale invariant contribution can be read directly from 
power counting in Eq. (js]) (leading to ^3 = C2)- Any dif¬ 
ferent scaling exponent can arise only from homogeneous 
solutions that do not need to balance the constant RHS. 
The scale invariant solutions of Eq. (||) live in a projective 
space whose dimension is lowered by unity compared to 
the most general form. These solutions do not depend on 
three separations but rather on two dimensionless vari¬ 
ables that are identified below. It will be demonstrated 
how boundary conditions arise in this space for which the 
operator B3 is neither positive nor self-adjoint. 

In Section 2 we present the transformation of variables 
that leads to a precise identification of the projective 
space. In this space we present the differential equation 
that needs to be studied, and derive the boundary con¬ 
ditions in the projective space. In section 3 we discuss 
the perturbation theory that leads to the solution of the 
scaling exponents of the zero modes in the limit (^h 0. 

It is shown that the choice of coordinates of section 2 
leads to a particularly transparent theory in this limit. 
In section 4 we present the perturbation theory in the 
limit 2 . It turns out that this is a singular per¬ 

turbation theory, and we discuss the analytic matchings 


across boundary layers and near the “fusion singularity” 
which are required to understand this limit, leading to a 
non-analytic dependence of ^3 on C,2- In Section 5 we deal 
with the non-perturbative calculation, culminating with 
solutions of C3 as a function of C,h throughout the range 
0 < C/i < 2 . It is demonstrated that the non-perturbative 
solutions are in agreement with the perturbative calcula¬ 
tions at the two ends of this interval. Section 6 is devoted 
to a summary and a discussion. 


II. TRANSFORMATION OF VARIABLES 


In this section we describe the transformation of vari¬ 
ables in the operator B3 to new variables that are denoted 
below as s, p, (j). We first note that equation (|^ is in¬ 
variant under space translation, under the action of the 
d dimensional rotation group SO((i), and under permuta¬ 
tions of the three coordinates. Accordingly, we may seek 
solutions in the scalar representation of SO(d), where the 
solution depends on the three separations ri2, r23 and r3i 
only. In the first stage we transform coordinates to the 
variables xi = |r2-r3p, a::2 = |r3-rip, X3 = |ri -r2p, 
defining 


T3(ri,r2,r3) = f3{xi,X2,X3) . (9) 


By the chain rule, 


9ijT3(ri,r2,r-3) = 2{ri3,d2f3 + ri 2 id 3 )f 3 {xi, X 2 , X 3 ), 

( 10 ) 


where du = (d/dr^J, 82 = {8/8x2), and ri2i = ru - r2i 
etc. 

Another application of the chain rule gives 

82 j 8 ii = 4ri3,r23j8i82 + iri2ir23jdi83 ( 11 ) 

—4ri3iri2j92i93 — 4ri2iri2ji93 — 2dij83 . 

(For brevity we display only the differential operators ex¬ 
plicitly). Using 2ri2 • ri3 = —xi -I- a;2 -I- X3, and similar 
identities, we can now obtain 


Sij 82 j 8 u = 2(a;i -I- a ;2 - a;3)i9i92 + 2(-a;i -I- a ;2 - X 3 ) 8 i 83 
+ 2{xi - X2 - X3)8283 - 4x38182 - 2d83 , (12) 


and 


ri2iri2j 


r 


2 

12 


82j8u 


-|-2(—a;i -|- X2 


{xi -X2- a;3)(a;i - a;2 + 3:3) 

- 8182 

a ;3 

X3)8i83 + 2 {xi - X2 - X3)8283 ( 13 ) 

4x38382 - 283 . 


Further calculations, using (0) and (|^) give 
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Bi2 = ri^2 


(d + Ch — ^)dij — Ch 


ri2iri2j 


d2jdu 


(14) 


„Ci ./2 


(d — I)oi + (d — I)(d + Ch)o2 + —03 

X3 


where 


oi = 2{xi + X2- X3)did2 + 2{-xi + X 2 - X3)did3 (15) 
+2{xi - X 2 - X3)d2d3 - 4X351 ^2 , 

02 = -253 , (16) 

03 = [a; 3 (a;i + X 2 - X 3 ) (17) 

-(xi - X 2 - X 3 )(xi - X 2 + a: 3 )] 5 i 52 . 


—p + cos{(j)) + -\/3siii((()), ^/ 3 cos{(j)) — sin((/)) 
—p + cos(^) — ■\/3sin((/)), —•\/3cos(^) + sin((/)) 
—/9 — 2 cos((/)), 2sin(0)/V3 


The Jacobian matrix J is substituted in the chain rule 
to give the transformation of the derivatives of /3 with 
respect to the x variables in terms of the new variables 
s, p, (j). It is convenient to perform the tedious calcula¬ 
tions using Mathematica ending up with expressions for 
oi, 02, and 03. 

We present the final long result as a table, in which 
each item is of the form 


The reader should note that the Oi operators do not de¬ 
pend on the parameters of the problem. The operators 
B23 and B31 are obtained from B12 by cyclic permuta¬ 
tions of the indices, thus completing the transformation 
of B3 to the X variables. 

Note that not every point in the xi,X2,X3 space cor¬ 
responds to a physical configuration. The triangle in¬ 
equalities between the pairwise distances translates to 
the condition 

2(xiX2 -I- X2X3 -I- X3X1) > xl + X2 + X3. (18) 

This inequality describes a circular cone in the xi, X2, X3 
space whose axis is the line xi = X2 = X3, tangent to the 
planes xi = 0, X2 = 0 and X3 = 0. The group of per¬ 
mutations between the Xi axes acts very simply on this 
cone, corresponding to a Cq operation. 

The presence of symmetries motivate a new parame¬ 
terization of the cone by three new coordinates s, p, (j): 

x„ = s[l — pcos((/)-I-|n7r)], (19) 

0 < s < 00, 0<p<l, 0 <^< 27 r. 


{s,p,(l>, n, TO, c} 
representing a term of the form 

cp- 

The trigonometric functions cos or sin appear if (/) is even 
or odd respectively. The Oi operators are sums of such 
terms. The tables listing the terms appearing in each of 
the operators are presented in Appendix A. 

The upshot of the transformation of the linear oper¬ 
ator B3 to the new coordinates is that we derive a sec¬ 
ond order linear partial differential operator in the s, p, (f) 
variables. At this point we take advantage of the scale in¬ 
variance of the differential equation which allows to seek 
scale invariant solutions of the form s‘’^^^f{p, 4 >). Act¬ 
ing on functions of this form, the operators sds and s^5g 
become scalar multiplications by and ^C3(5C3~1) re¬ 
spectively. The action of the operator B3 yields an equa¬ 
tion for /(p, (j)), which is the basic equation we study in 
this paper 


The new space is a direct product of three intervals, and 
this fact will simplify the discussion of the boundary con¬ 
ditions. The s coordinate measures the overall scale of 
the triangle defined by the original coordinates, and 
configurations of constant p and (j) correspond to simi¬ 
lar triangles. The p coordinate describes the deviation 
of the triangle from the equilateral configuration [p = 0) 
up to the physical limit of three collinear points attained 
when p = 1; (j) does not have a simple geometric mean¬ 
ing. Finally we note that the variables s,p and cos(3(/)) 
are symmetric in the Xi variables. Accordingly any func¬ 
tion of these variable is automatically invariant under the 
permutation of the Xi variables. We will use this property 
below. 

The final form of the equation is achieved by trans¬ 
forming Oi operators to the variables s,p, 4 >. To this end 
we compute the Jacobian of the transformation (H) us¬ 
ing Mathematica, 

a(s, 1 

5(Xi,X 2,X3) V3s2p 


^3(C3)/(p, <)>) = [a{p, 4>)dl + b{p, (j))dl -f c(p, 4>)dpd^ (21) 
+u(p, 4 >, C3)5p -f v{p, (j), (3)8^ + w{p, (j), C3)] f{p, (t>) = 0 ■ 

We note that Eq. iE) can be written in a coordinate-free 
form as 

[-V-P(p,</.)-V + q(p,</.,C3)-V (22) 

+ w(p,^!),C 3 )]/(p,^) = 0 , 

where V is the gradient operator in the p, (j) space. The 
identification of the tensor P and the row vector q is ob¬ 
tained by comparison with the explicit form (pl|). The 
new operator B3 depends on C.3 as a parameter and it acts 
on the unit circle described by the polar p, (j) coordinates. 
The circle represents the projective space of the physical 
cone described above. We will see that the availability of 
a compact domain (the projective space) will lead to the 
existence of a discrete spectrum of the zero modes. 

The discrete permutation symmetry of the original 
Eq. (^ results in a symmetry of Eq. with respect 
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to the 6 element group generated by the transformation 

^ + 27 r /3 (cyclic permutation of the coordinates in 

the physical space) and ^ (exchange of coordi¬ 
nates). This symmetry extends to a full 17 ( 1 ) symmetry 
in the two marginal cases of C/t = 0 and C/t = 2 (see 
for a discussion of the latter limit) for which all the coeffi¬ 
cients in ( pl| ) become 0 -independent. The coefficients in 
iE) all have a similar structure, and for example a(p, 0 ) 
reads 

a(p, 0 ) = ^[1 - pcos (0 -I- | 7 rn)](‘”'*“^)/^d(p, 0 -I- | 7 rn) , 

n 

where d(p, 0 ) is a low order polynomial in p, cos 0 and 
sin 0 which vanishes at p = 1,0 = 0 . We see that the 
coefficients are analytic everywhere on the circle except 
at the three points p = 1 , 0 = 27 rn /3 where n = 0 , 1 , 2 . 
These points correspond to the fusion of one pair of co¬ 
ordinates, and the coefficients exhibit a branch point sin¬ 
gularity there. This singularity is a consequence of the 
non-analyticity of the driving velocity field whose eddy- 
diffusivity is therefore characterized by a non-integer ex¬ 
ponent. This singularity leads to a nontrivial asymptotic 
behavior of the solutions which had been described before 
in terms of the fusion rules l9|Jl^. Since the coefficient 
d vanishes at the fusion point, the singularity is always 
sub-leading with respect to terms that come from non¬ 
fusing coordinates. Indeed, in [pUl^ it was explained that 
exposing the singularity calls for taking a derivative with 
respect to the fusing coordinates. Note that for 0/1 = 2 
the singularity disappears trivially. For Ch = 0 there is 
also no singularity since d exactly compensates for the 
inverse power. 

The boundary conditions follow naturally when one 
realizes that is elliptic for points strictly inside the 
physical circle. In the presentation of Eq. ( p^ this means 
that 

det(P) >0 for p < 1 . ( 23 ) 

This property is a consequence of the ellipticity of the 
original operator B3. On the other hand B3 becomes 
singular on the boundary p = 1 , where the coefficients 
a(p, 0 ) and c(p, 0 ) vanish. In other words, 

P-n|p=i= 0 , ( 24 ) 

where n is a unit vector normal to the boundary. This 
singularity reflects the fact that this is the bound¬ 
ary of the physical region. It follows that B3 re¬ 
stricted to the boundary becomes a relation between 
the function /(p=I, 0 ) = p( 0 ) and its normal derivative 
5 p/(p, 0 )|p=i = h{cj)). The relation is 

4 ’) + uh{l, 0 ) -t vg'{l, 0 ) -b wg{l, 0 ) = 0 . ( 25 ) 

Solutions of Eq. ( pT[ ) which do not satisfy this bound¬ 
ary condition are singular, with infinite p derivatives at 


p = 1 . Such solutions are not physical since they in¬ 
volve infinite correlations between the dissipation (sec¬ 
ond derivative of the field) and the field itself when the 
geometry becomes collinear, but without fusion. 

It is important to stress that the reason for the regu¬ 
larity of the solutions that do satisfy the boundary con¬ 
ditions is that detP vanishes as a simple zero near the 
boundary, i.e. like (I — p). Consequently it is possible to 
find solutions that behave like (I — p)° ^ 0 ( 1 ) near the 
boundary. 

Finally, we are facing the problem of solving Eq. © 
together with the boundary condition (p^). This homo¬ 
geneous boundary value problem always have a trivial 
solution / = 0 . The condition for the existence of a non¬ 
trivial solution is that P3 is not invertible. We expect 
that non-invertibility will not be a generic phenomenon 
for an arbitrary value of 03 . Einding the set of 03 for 
which P3 is not invertible becomes a generalized eigen¬ 
value problem with the exponents 03 playing the role of 
eigenvalues. We first discuss this program in the limit of 
small 0/j. 

III. PERTURBATION THEORY NEAR 0^ = 0 

The shape of our compact domain invites a Fourier 
representation in 0 for the function /(p, 0 ). The per¬ 
mutation symmetry implies that only cos will appear in 
this representation, and the index will be divisible by 3 . 
The general equation ( pl| ) will then mix different Fourier 
modes. On the other hand, in the limit 0 /i = 0 the situ¬ 
ation simplifies considerably. 

A. Zero Modes at 0/, = 0 

The differential equation reduces in the case 0 /i = 0 to 
the simple form 

[p 2 (l - p^)dl + p(l - dp^)dp + dl + Ap^]/(p, 0 ) = 0 , 

( 26 ) 

where 

\=\U\ 4 z + d-l) . ( 27 ) 

Since the coefficients of Eq. (|^ are independent of 0 
the Fourier modes are decoupled, and we may seek so¬ 
lutions of the form /m(p)cosTO 0 with m divisible by 3 . 
The functions fm obey the ODEs 

[p^(l — p^) 5 p -f p(l — dp^)dp — m? + Ap^]/m(p) = 0 . 

( 28 ) 

Examining the resulting differential equation we note 
that the coefficient of the highest derivative in p van¬ 
ishes as a double zero at p = 0 and as a single zero at 
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p = I. Since the order of the zero is not more than the or¬ 
der of the derivative the Frobenius theory |^l| of regular 
singularities is applicable to both boundaries. Namely, 
the complete family of solutions in the vicinity of a given 
singular point po is spanned by functions that can be 
represented as 

CXD 

/m (p) = {p- PoY' log{p - Pof' - poY, 

P—0 

i = l, 2 , ki = 0 orl, (29) 


where the sum is convergent in a neighborhood of po- 
One of the indices ki is always zero. When the indices 
zi,2 are different, we chose arbitrarily zi > zi and then 
ki is zero. When jzi — Z2I is not an integer one also has 
^2 = 0 . In cases in which the indices Zi coincide we will 
chose fci = 0 , and ^2 = 1 . The numerical values of the 
indices are obtained by substituting a solution (p — po)^ 
in the differential equation, collecting the coefficients of 
the leading terms near the singularity, and equating to 
zero. We refer to the indices Zi as the “Frobenius expo¬ 
nents” . Eq. (^ ) has regular singularities at both bound¬ 
aries p = 0 , 1 , with Frobenius exponent sets of —m, m 
and 0 , {3 — d )/2 respectively. The singularity at p = 1 
arises from the singularity of the original PDE ( ^l|) at this 
boundary, and the boundary condition picks the regular 
solution, i.e., fm{p) (1 ~ p)° as p —> 1 . On the other 
hand, the singularity at p = 0 is an artifact of the trans¬ 
formation to polar coordinates; however, analyticity of 
the solution at p = 0 requires that fm{p) ~ p™ as p ^ 0, 
specifying the second boundary condition for ( 28 ). 

To proceed with the solution of Eq. (|^) we make the 
transformations fm{p) = p"^fm{p) and p = p'^, and ob¬ 
tain the hypergeometric equation for fm, 

{(1 - p)p 5 | + [m + 1 - {m + + l)]p)dp ( 30 ) 

-h|[A - m{m - 1 -h d)]]fm{p) = 0 . 

This equation is standard, see |^. The unique solution 
of Eq. (| 30 | ) which satisfies the boundary condition at 0 
(up to a multiplication by a constant) is 

fm{p) = 2Fi{a,b;m+l-,p), ( 31 ) 

where 

a=i(TO-iC 3 ), b = ^{m +d -1 + ^Cs) ■ ( 32 ) 

It follows from the theory of the hypergeometric func¬ 
tions that the function defined in (^) is regular at 
p = 1 only if either a or b equal —n, where n is a nonneg¬ 
ative integer. In such a case fm becomes a polynomial of 
degree n. The spectrum of Ca now follows from (^) and 
consists of two sets, 

CaV.n) = ^Krn,n) = - 1 + m + 2n). 


(33) 

Since A(C3+(„_„)) = A(C3-(„_„)), [cf. Eq. (|^], the corre¬ 
sponding eigenfunctions depend only on n and m. Ex¬ 
pressed in terms of the original variables the eigenfunc¬ 
tions are 


f(m,n){pA) (34) 

= p™ 2Fi (—n, n-|-TO-|- {d— l)/2;m-|- l;p^) cos(m0) . 

The zero-modes of the first set, with positive values 
of Ca are the ones to be matched at the outer scale. Of 
these, the most relevant non-trivial zero mode is /o,i with 
1) = 4 , still less relevant than the scaling of the forced 
solution. 


B. Zero Modes for 0 < ^ 1 


Perturbation theory for Ch 0 will be carried out by 
expanding the solution for small Ch in terms of the ((h = 0 
eigenfunctions in a procedure very similar to time 

independent perturbation theory in quantum mechanics. 
We first cast Eq. (^^ in Sturm-Liouville form 


.p(l_p 2 )(d-a )/2 


do 


p(i-P^)("-i)/^a, 


(35) 


“ = -^fmip) ■ 

It follows that the zero-modes f(rn,n) form an orthogonal 
set with respect to the inner product 


if, 9)0= J -^PiX- P'^)^'^ ^^^^f{P^ 4 ') 9 {pA) ■ (36) 

p<i 


We now assume that we may expand the eigenfunction 
and eigenvalues as 

f{p, </>) = hm,n){p. </-) + (P, 4>) + 0(CD , (37) 

Ca = (,3{m,n) + C 3 ^ + O(Cft) ; (38) 

writing 


FiChXs) — B{0, Ca(m,r!,)) + Ch dchF{0, Ca(m,n)) 


(39) 


+ C 3 C3(m,n)) + 0{Ch) 


Eq ( 21 ) becomes, to order (h, 


F{0, C3im,n))f^^\p, 4>) + d^f^B{0, C3(m,n)) 


(40) 


+C 3 ^d^^B{0X3{m,n)) /(m,n)(Pi <(>) — 0 ■ 


Taking the inner product ()q of Eq. ( ^ ) with /(m,n)j and 
using the self-adjointness of B{ 0 , C3(m,n)) with respect to 
this inner product, there follows an expression for the 
first correction to C3, 



m,n m,n )) Q 
m,n ))/(m,n))Q 


(41) 


The integrals in (© were programmed using Mathemat¬ 
ical and we used the program to generate expressions for 
the first few values for ^3^^ presented below: 
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TO 

n 

*’3(m, 

0 

0 

0 

0 

1 

4 

0 

2 

8 

0 

3 

12 

3 

0 

6 

3 

1 

10 

6 

0 

12 


p+(i) 

S3 

^3(m, 


^-(1) 

S3 

.n) 

0 

-2d+2 

-1 

2(2-£i) 

d-1 

-2d 

- 2 

3-d 

d-1 

2 (159-1-91 d-31 d^-19 d^-2 d*) 

-2d 

- 6 

-3 (90+62 d-10d^-9d^-d*) 

-48+4 d+32 d2_|_ii d3+d4 

-48+4 d+32 d2 + ll ds+d* 

6 (780+503 d-89 d^-103 

-2d 

- 10 

4200+3010 d-210 d^-476 d^-93 d'‘-5 d® 

-480-8 d+324 d^ + 14:2 d3+21 d*+d^ 

-480-8 d+324 d2 + 142 d3+21 

-3 (15+8 d+d^) 

-2d 

A 

(5+d)(7+2d) 

-2+d+d^ 


2-d-d^ 

(9+d) (-489-131 d+136d^+53 d^+5d^) 

-2d 


(9+d) (441 + 135 d-104d^-42 d^-Ad"^) 

-144-36 d+100 d^+65 d^ + 14 d^+d^ 


-144-36 d+100 d^+65 d-^ + lA d^+d^ 

3 (3+d) (11+d)^ (37O+2I6 d+37 d^+2 d^) 

-2d 

- 10 

5 (11+d) (6558+6508 d+2409 d^+412 d^+33 d^+d^) 

3840+544 d-2584 d'-‘ - 1460 d^ -310 d'* -29 d^ -d« 

-3840-544 d+2584 d^ + 1460 d^’+OlO d*+29 d^+d" 


(42) 


IV. PERTURBATION THEORY IN THE LIMIT 

Ch^2 

In this section we treat the case C2 = 2 — 1 per- 

turbatively. The perturbation theory is singular, since 
the leading order approximation, obtained by setting the 
small parameter ((2 = 0 is not valid throughout the en¬ 
tire domain. We employ boundary layer techniques and 
matching of asymptotics to obtain an asymptotic approx¬ 
imation for ((3 in this limit. As expected, ^3 goes to 0 with 
C2) but with a non-trivial dependence 

Ca = 0(\/C2/logC2) • (43) 


A. Leading order solution 


As in the case = 0 discussed above, substituting 
C/i = 2 into the zero-mode equation (H) yields an equa¬ 
tion with coefficients which are independent of the vari¬ 
able (j). This equation reads 


[p 2 (i _ p2^2g2 ^ _ ^ 2 )(i _ ^^2 ^ ^^^ 2 )^^ ( 44 ) 

+ (1 - p'^)dl + wip)]fiP, ^) = 0, 

Mp) = P^Yd + 2)(d - 1) + (| - i)(('^ - i)(p" + 1) 

- l)l . (45) 


that the Frobenius exponents at 0 are ±m, and choose 
the solution which behave as p™ at 0, thus providing the 
boundary condition at 0. 

The behavior near the boundary p = 1 is different. 
In contrast with the case Q < C,h < 2 where one of the 
Frobenius exponents is always 0, when Ch = 2 the Frobe¬ 
nius exponents depend on the value of Ca, and generally 
neither is 0. Thus, the qualitative behavior of the zero 
modes for 0 < ((2 = 2 — <C 1 near p = 1 is differ¬ 
ent from that of the solutions of Eq. (^) where C,h is set 
equal to 2. In other words, we expect the existence of a 
boundary layer near p = 1; Eq. (|^) describes well only 
the behavior of the outer solution, namely the leading 
order approximation away from the boundary layer. 

The outer solution may now be written explicitly. 
Using standard transformations Eq. ( ^ ) may be 
reduced to a hypergeometric equation, whose solution 
which obeys the p = 0 boundary condition is 

fmip) = P™ 2E1 (to/2, (to + 1)/2;to + l;p^) . (48) 

Since we anticipate that ^3 <C 1, we give in ( ^^ the solu¬ 
tion of equation (|^) taking ^'3 = 0. This approximation 
is justified to leading order a-posteriori, when we find 
that /'3 is indeed small. The outer solution is constructed 
from the /m’s by 

/°“‘(p,(/>)= ^ Vmfmip) cos mcj) (49) 

m—0,3,6,... 


This significant simplification is a consequence of the 
higher symmetry of the passive scalar equation in this 
limit, as discussed in detail in Ref . |0,||. 

Making use of the symmetry, we look for solutions of 
the form 

f {p A) = fm{p) cos m4> TO = 3n, n = 0,1,2 ... (46) 

The functions fm satisfy the ODEs 

[ p 2 (l _ ^ 2 ) 2^2 ^ ^(4 _ ^ 2)(4 _ 2^2 ^ 

p^)+w{p)]fm{p)=0 . (47) 

Equation ( ^^ has regular singularities at the points 0 
and 1. In close analogy with the case C/, = 0 we find 


The coefficient Vm are unknown at this stage and they 
will be determined by matching with the inner solution. 

The outer solution is a valid asymptotic approxima¬ 
tion of the true solution only when C2 1 ~ P- Since we 
expect the solution to vary very rapidly within a bound¬ 
ary layer near p = 1, the inner solution, valid within the 
boundary layer, should be expressed in terms of a “fast” 
variable r that changes on a scale inversely proportional 
to C2, namely, 

p=l + C 2 T. (50) 

We therefore change the variable p to r in the differen¬ 
tial equation and keep the leading terms in (^ 2 - The 
resulting equation for the inner solution is 
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{t[t- 


d + 1 


■pi^M + ( 


(d + 3)(d-l) 
4d 


PW - -jT)dr 


+w(l)/4}/“(r,(/)) = 0 , 


(51) 


where 


where 


/—pr (d+3)(<i-l) l 

^ [ 2{d+i) J d — 1 

■n r(d+3)(^*—1) I 11 4 

^ L 2(d+i)—^ 2 J 



(57) 


Pi^) = ^ log(l — cos(/)')(cos(()'— cos2(()'), (52) 

0±27r/3 

and the function w is defined in (^5|) . It is crucial for the 
matching between inner and outer solution to realize that 
Eq. (H is not valid near the fusion points (j) = ^nrijS. 
The significance of this fact is explained below. 

Equation & has a regular singularity at r = 0 which 
corresponds to the p = 1 boundary, and as in the other 
cases demanding regularity gives a boundary condition 
at this point. The equation can be again transformed 
into the hypergeometric equations, so it has solutions of 
the form 

/“('r,</>) = P(</') (53) 

p ( jd-m l, (d + 3)(d-l) 2d _ 

^ ^ ^ V 4 ’2’ 2(d+l) ’ d+lp(<())y 

Here again we used the smallness of (3 and neglected 
terms of higher order in it. The function f™ provides an 
asymptotic approximation to the actual solution when 
—T (1/C2)- It depends on the function ^(0) that will 
be determined by matching to the outer solution. 


B. Asymptotic matching 

The next step involves matching of the two approxi¬ 
mations /™ and in their common region of validity 

C2 < 1 - p < 1 . (54) 

We perform the matching using standard boundary layer 
methods by examining the asymptotic behavior of f™ 
and /°“* in the matching region ( p^ and balancing co¬ 
efficients. We make use of the asymptotic behavior of 
the hypergeometric functions . For the outer solution 
we need the asymptotics of the hypergeometric function 
with argument close to 1, yielding 


\p,^) - 

i-P<i ™ 


1— m\/ 2 {l — p) cos(m0) . 


(55) 


For the inner solution we use the behavior of the hyper¬ 
geometric function for large negative values, and obtain, 
using the relation between p and r 


r(p,</>) 


- p((/>) 


i-p<i 


1 + cCs 



(56) 


To proceed with the matching of ( ^ and (^^, we first 
match the 0(1) terms [for small (1 — p)], giving 

piij)) = E cos(m^), (58) 

m 

i.e. T^Vrn are just the coefficients of the Fourier series 
of Next we have to match the 0{^/l — p) terms 

giving 


cCa 


--pW = E —mV2 2'^Vm cos{m(j)) . (59) 


^/C,2p{4>) 

We expand p(^)“^/^ in Fourier series, 
p(<^)-i/2 = cos{m(j)) 


(60) 


and substitute in ( b9D to get, using 

cCa 


E 7 r^PnVm{cos[(n - to)((>] -b cos[(n -b m)^]) ( 61 ) 
— ^ zy C2 

mn ^ ’ 

= E cos{m(j)) . 


Equating Fourier coefficients of the same order yields fi¬ 
nally 


E 


cCs 


Pn{^m+n “b —n|) — mi^m- 


(62) 


The matching condition (^) is a generalized eigenprob- 
lem for the infinite vector Vm, with eigenvalues ^3. 

Since ^3 appears in (|^) only through the combination 
Ca/v^i one would be led to conclude that ^3 = O(v^) 
for all non-trivial solutions of ( |6^ with different numeri¬ 
cal coefficients. This consideration is modified, however, 
since the coefficients Pn themselves also depend on (2 as 
will be now demonstrated. 

It follows from the definition of p{(j)) (|^, that 


p{(l)) = O[0^ log(0)] (63) 


for (j) small. Note that the Fourier coefficients p„ are 
written as integrals 


i 

^ h \/p{4>) 


for n > 1 . 


(64) 


All these integrals diverge at 0 = 0, which is precisely 
the fusion point where the approximation (53) ceases to 
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be valid, requiring a more careful examination of the be¬ 
havior in this region. In deriving Eq. we made the 
approximation 

(1 — p cos ~ 1 — 5C2 log(l — p cos (j)) 

- 1 - iC2log(l - cos(/i) , (65) 

for 1 — p = 0 {(2)- When (j)^ = 0{C,2) the second approx¬ 
imation is no longer valid, and one has instead 

(1 - pcos(/))“‘^='/^ - 1 - iC2log[(l - p) + , (66) 

so that instead of ( [^ ) we get a similar equation in which 
is replaced by log(—C2T + The resulting equa¬ 
tion is no longer integrable in terms of hypergeometric 
functions. However, we still expect the asymptotic ap¬ 
proximation ( ^ ) to be valid, but for small (j) the function 
p((^) is no longer given by (^. We may estimate p{(j)) for 
small (j) by the following consideration. Examining Eq. 
(H we see that p{(j)) is the value of r where the coeffi¬ 
cient of the second derivative crosses over from quadratic 
behavior to linear behavior, and the coefficient of the first 
derivative crosses over from linear to constant behavior. 
When (j)^ = 0(^2) there is also a crossover, but to a lin¬ 
ear or constant function times a logarithmic function of r. 
The logarithmic function changes very slowly, and may 
be approximated roughly by a constant. The crossover 
occurs when 

|r| ~ C 2 ^ 0 g{C 2 \T\ + (j}^) ~ C2log(C2) ■ (67) 

We are thus led to make the following approximation 

p(^) - C2logC2, for (j)'^= 0 {C 2 )- (68) 

This estimate implies that when calculating the Fourier 
coefficients in (^) the integral should be cut off at (/) ~ 
giving the coefficients a |fog(C2) dependence on 
((2 • It is now possible to balance powers in the matching 
Eq. (p^, obtaining the order of magnitude relation (^). 


V. SOLUTIONS FOR GENERAL VALUES OF (h 


for discrete values of Ca for any given values of Cft, and the 
dimensionality d. One solution is known to exist always, 
a constant zero-mode associated with (C3 = 0. Our aim 
is to find the lowest lying positive real solutions Ca for 
which the determinant vanishes. 


The discretization of the operator B 3 was carried out as 
follows: We defined a nine-point finite difference scheme 
for the second order equation (|^). The discretization of 
the boundary conditions at p = 1 (^) is achieved us¬ 
ing the same scheme, which requires in this case only 
three boundary points and one interior point, since on 
the boundary the radial derivative appears in first or¬ 
der. Using the symmetry of the problem we restricted 
the domain to one sixth of the circle, 0 < (j) < tt/S. The 
symmetry implies that original problem on the full circle 
is equivalent to the problem on the reduced domain with 
simple Neuman boundary conditions ^^/(p, 4>) = 0. on 
the new boundary lines 0 = 0, tt/S. As explained above, 
the discretized problem is a matrix eigenvalue problem 
HafF = 0, where B 3 is a large sparse matrix, whose rank 
depends on the mesh of the discretization, and ^ is the 
discretized version of the zero-mode /. We used NAG’s 
sparse Gaussian elimination routines to find the zeros of 
det(i?3), and determined the values of Cs for these zeros 
as a function of (h ■ The results of this procedure for space 
dimensions d = 2, 3,4 are presented in Figs. 1,2, and 3. 
The zero modes that correspond to any given value of 03 
can be found straightforwardly by inverse iterations of 
the matrix B 3 with an arbitrary initial vector. 


A. Results 


The various branches shown in Figs. 1-3 can be orga¬ 
nized on the basis of the perturbation theory near C/i = 0 
which was presented in Section 2. At Cft, = 0 we identify 


For general values of and d the differential equa¬ 
tion (il , which has variable coefficients is not accessible 
to analytic techniques. In this section we present numer¬ 
ical solutions of the scaling exponents ((3 for arbitrary 
values of C/i, using a discretized version of the operator 
B^. Since the differential problem is a linear homoge¬ 
neous equation with linear homogeneous boundary con¬ 
ditions, the discretized problem is also a homogeneous 
linear equation, implying that non-trivial solutions exist 
only when the determinant of the discretized operator 
vanishes. This determinant depends parametrically on 
03. Since the differential operator is defined on a com¬ 
pact domain we expect the determinant to vanish only 

FIG. 1. The scaling exponent 03 as a functions of C,h found 
as the loci of zeros of the determinant of the matrix , for d=2. 


FIG. 2. Same as Fig.l, but for d = 3. 

FIG. 3. Same as Fig.l, but for d = 4. 

the actual starting points of the branches with the ana¬ 
lytic solutions for the lowest lying positive values of 03, 
which are 4, 6,8 etc. In addition, for d = 2 we observe the 
highest negative value which is —2. Measuring the slopes 
of the various branches at 0?i = 0 we find full agreement 
with the perturbative predictions. 

We see that in all dimensions the branch which begins 
at 0?i = 0, 03 = 4 continues in the non-perturbative re¬ 
gion without crossing any other branch until it ends at 
0^ = 2, 03 = 0. This branch is a continuation of the 
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lowest lying positive branch predicted by the perturba¬ 
tion theory. The negative branch (shown only for d = 2 ) 
never rises above its perturbative limit and is not rele¬ 
vant for the scaling behavior at any value of Ch ■ Note also 
that the point = 2, ^3 = 0 is an accumulation point 
of many branches, and we display only a part of the ac¬ 
tual spectrum near this point. These branches seem to 
approach the accumulation point with a slope that grows 
without limit. This finding is in agreement with the an¬ 
alytic result of the perturbation calculation presented in 
Section 4 . 

VI. CONCLUSIONS 

It is well known by now that there exists a disagree¬ 
ment between the scaling exponents ^4 and the higher or¬ 
der exponents C.n computed via perturbative approaches 
and the predictions of another approach based on the 
fully fused structure functions. The latter approach 
seems to be consistent with the results of numerical sim¬ 
ulations in two spatial dimensions. The main conclusion 
of the present paper is that this disagreement cannot be 
ascribed to a formal failure of the perturbation theory. 
If we accept the statement that the scalar diffusivity is 
irrelevant, and compare the predictions of perturbation 
theory at both ends of the range of the allowed values of 
the parameter (^h, we find that they are in excellent agree¬ 
ment with the non-perturbative calculation of the scaling 
exponent, again subject to the assumption that the diffu¬ 
sivity is irrelevant. Therefore, if we want to understand 
the discrepancy between the two approaches mentioned 
above, there are a few possibilities that have to be sorted 
out by further research: 

(i) The crucial assumption that goes to the fully fused 
approach, which is the linearity of the conditional aver¬ 
age of the Laplacian of the scalar, is wrong. 

(ii) The computation of the zero modes which is achieved 
by discarding the viscous terms in Bn is irrelevant for 
the physical solution. It is not impossible that the limits 
C/i —*■ 0 and K ^ 0 do not commute, giving rise to some 
wicked properties of the very small (^h regime. That this 
is a possibility is underlined by recent calculations of a 
shell model of the Kraichnan model , in which it was 
shown that the addition of any minute diffusivity changes 
the nature of the zero modes qualitatively. 

(iii) Lastly, and maybe most interestingly, it is possible 
that the physical solution is not strictly scale invariant 
through all the range of allowed distances j^. In other 
words, it is possible that lA3(ri,r2,ra) is not a homo¬ 
geneous functions with a fixed homogeneity exponent ^3, 
but rather (for example), that ^3 depends on the ratios of 
the separations (or, in other words, the geometry of the 
triangle defined by the coordinates). If this were also the 
case for even correlation functions J-2n-, this would open 
an exciting route for further research to understand how 


non-scale invariant correlation functions turn, upon fu¬ 
sion, to scale invariant structure functions. 

In light of the numerical results of Ref. and the ex¬ 
perimental results displayed in [|l 6 [|l 7 | ] we tend to doubt 
option (i). If we were to guess at this point we would 
opt for possibility (ii). More work however is needed to 
clarify this important issue beyond doubt. 
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